Gate-controlled conductance through bilayer graphene ribbons 
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We study the conductance of a biased bilayer graphene flake with monolayer nanoribbon contacts. 
We find that the transmission through the bilayer ribbon strongly depends on the applied bias 
between the two layers and on the relative position of the monolayer contacts. Besides the opening 
of an energy gap on the bilayer, the bias allows to tune the electronic density on the bilayer flake, 
making possible the control of the electronic transmission by an external parameter. 
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The prospective use of graphene in nanoelectronics re- 
quires the possibility to open gaps in its band structure 
in a controllable way. Due to the chiral nature of the car- 
riers £Q , it is not easy to open gaps and confine carriers in 
a single graphene monolayer. Carbon-based structures as 
nanoribbons [2| [3] , nanotubes [IHS] , and graphene bilay- 
ers [7HT2"] are viable materials for nanoelectronics, since it 
is feasible to change their electronic characteristics from 
semiconducting to metallic as a function of geometric or 
external parameters. Bilayer graphene is a good candi- 
date because a gap can be opened and controlled by an 
applied bias between its two layers. Monolayer graphene 
nanoribbons (MGNs) stand out as optimal electrodes 
for systems based on bilayer graphene, with the aim of 
achieving the best integration of nanoelectronic compo- 
nents. Thus, it is important to study the electronic trans- 
port of bilayer graphene nanoribbons (BGNs) with MGN 
contacts. Previous work has focused on the electronic 
transport through bilayer graphene flakes in absence of 
external gates [13] ■ In such a case the conductance shows 
strong oscillations as a function of the energy of the in- 
cident electron and the length of the bilayer region. In 
this work we show that the conductance of BGNs con- 
nected to MGNs strongly depends on the way the bilayer 
is contacted and on the applied gate voltage. This allows 
for an external control of the electronic properties of the 
system. 

Geometry. We analyze electronic transport in the lin- 
ear regime through a gated BGN connected to two metal- 
lic MGN contacts. The monolayer leads can be either 
armchair or zigzag graphene nanoribbons [3] |3J, serv- 
ing as contacts to armchair or zigzag bilayer flakes re- 
spectively. In both cases two configurations are possible: 
the bottom-bottom (1— sT) and the bottom-top (1— >• 2), 
where the ribbon leads are connected to the same (1—^1) 
or to a different monolayer (1— > 2) of the bilayer flake. 
We consider BGN of width W and length L, and restrict 
our study to narrow nanoribbons in the energy range 
for which only one incident electron channel is active. 
The bias is applied symmetrically with respect to the 
top (-V/2) and the bottom (V/2) layers. 

Electronic structure of constituents. The band struc- 
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FIG. 1. (Color online) Schematic bandstructure of a biased 
graphene bilayer with and without interlayer coupling 71 . (a) 
Biased bands without interlayer hopping. A positive (nega- 
tive) bias V/2 is applied to the bottom (top) layer, producing 
a rigid shift of the corresponding linear dispersion relation 
(red dotted lines for the bottom layer, blue dashed lines for 
the top one), (b) When 71 is switched on, gaps open at the 
band crossings, yielding the well-known Mexican hat shape 
and split-off bands of the bulk bilayer bandstructure. This 
picture allows to identify the top (T)/bottom (B) character 
of the different branches of the bilayer dispersion relation. 

ture of graphene has two inequivalent valleys. Within 
one valley, the low energy properties of graphene are 
well described by the two dimensional Dirac equation, 
H = VF&-p, where vp ~ 1 x 10 6 m/s is the Fermi velocity, 
p is the momentum operator relative to the Dirac point 
and o~i are the Pauli matrices. The Dirac Hamiltonian 
acts on a two-component spinor, (4>a,4>b)i representing 
the amplitude of the wavefunction on the two inequiv- 
alent triangular sublattices of graphene, labeled A and 
B. The band structure of armchair graphene nanorib- 
bons is obtained from the Dirac equation with the ap- 
propriate boundary conditions 3 . In all nanoribbons 
the transverse momentum is quantized. For the armchair 
MGN case, when the number of carbon atoms across the 
width of the ribbon is equal to 3to + 2, being m a pos- 
itive integer, the smallest transverse momentum is zero. 
This yields the ribbon metallic with an energy dispersion 
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vfPx, where p x is the momentum along the nanoribbon. 
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FIG. 2. (Color online) Schematic band structure of the bilayer 
armchair nanoribbon with monolayer contacts. Upper and 
lower panels show the bottom-bottom (1— >1) and bottom-top 
(1— >2) configurations. Left and right dispersions correspond 
to the metallic monolayer ribbons acting as electrodes, and 
the central dispersion corresponds to the biased bilayer rib- 
bon. The applied gate bias between bottom and top layers 
is V. For widths for which armchair monolayer ribbons are 
metallic, the band structure of the bilayer is not affected by 
the confinement. In the band structure of the bilayer we in- 
dicate the three relevant gaps, V, Ai, and A2. 



In addition to confined states, zigzag MGNs support 
zero energy surface states located at the edges of the 
ribbon 2 . In reciprocal space, surface states occur be- 
tween the two Dirac cones and their number cannot 
be described by the two dimensional Dirac equation, 
which is only valid for the low energy physics near the 
cones. Therefore, to describe zigzag graphene nanorib- 
bons we use a nearest-neighbor tight-binding Hamilto- 
nian H = — t^2(afbj + h.c). Here di{bi) annihilates an 
electron on site i of sublattice A(B), and the hopping pa- 
rameter t is related to the Fermi velocity by vf — g a *> 
where a is the graphene lattice constant, a ~ 2.46 A. 

Bilayer graphene consists of two coupled graphene lay- 
ers with inequivalent sites A\,B\ and A2, B2 on the bot- 
tom and top layers respectively. We consider the Bernal 
stacking in which the B2 sublattice is exactly on top of 
the sublattice A\. Within one valley, the low energy 
properties of a biased bilayer graphene are well described 
[T4] by the Hamiltonian 

V 71 

H B G = VFTo®5-p+ — T z ®(T + — {T x ®(T x -Ty®(Ty) : (1) 

where 71 ~ t/10 is the hopping parameter between the 
closest carbon atoms belonging to different layers, <7j 
are again the Pauli matrices for the sublattice degree 
of freedom and r» are the Pauli matrices for the layer 
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FIG. 3. (Color online) Same as Figj2] but for zigzag termi- 
nated nanoribbons. For zigzag monolayer nanoribbons, states 
in different valleys transport charge in opposite directions. In 
the biased bilayer ribbon, central part, edge localized states 
with energies ±V/2 appear connecting the Dirac points. At 
the Dirac points of the bilayer ribbons, topological surface 
states exist in the gap. In wide ribbons these surface states 
close the energy gap. In narrow bilayer ribbons states in dif- 
ferent edges interact and open a gap A in the spectrum. 



index (oo and To are the unit matrices in both sub- 
spaces). This Hamiltonian acts on the four-component 
spinor {4>\, 4> B ,ip\ 1 4>' B ) representing the amplitude of the 
wavefunction on the two sublattices A and B of the 
two layers 1 and 2. The energy bands are £%q(p) = 
2 2 1 v 2 



t + \ ± |^4v|.p 2 (V 2 +7?)+7j. The low en- 
ergy band has a Mexican hat shape with a minimum gap 
Ai = 7i|V|/(\/7i + V 2 ), see Fig. |l| The mi nimum gap 

of the second subband is A2 = 2\/7 2 + O^) 2 and occurs 
at p = 0. Fig. [T] illustrates how the Mexican hat shape 
and the split-off bands arise: without interlayer hopping, 
the applied bias shifts the linear band dispersions of the 
two layers; the interaction between layers opens gaps at 
the intersections of the bands. 

The electronic structure of an armchair BGN depends 
on the width of the ribbon. As in the monolayer case, 
when the number of carbon atoms along a BGN layer 
is equal to 3m + 2, the smallest transverse momentum is 
zero, and the dispersion of the armchair BGN is 6bg(Px), 
p x being the momentum along the ribbon, see Fig. [2] 
In the case of zigzag biased BGNs the system supports 
two kinds of surface states [TS]: (i) states with energies 
~ ±V/2, similar to those occurring in zigzag MGNs, and 
(ii) valley-polarized states with energies in the gap. At 
each edge of the ribbon there are two surface states car- 
rying current in opposite directions and belonging to dif- 
ferent valleys. These states have a topological nature [TB] 
but the metallicity of the edge is not protected against 
intervalley scattering nor against interedge intravalley 



scattering. The latter occurs when the ribbon width is 
smaller than the penetration length of the surface states 
[17], I ~ y/3^a, that for realistic values of the inter- 
layer hopping is around 17a. This large £ value produces 
an interedge scattering gap A in the spectrum of narrow 
zigzag BGNs, see Fig. [3J Although the valley-polarized 
surface states can be modeled with the Dirac Hamilto- 
nian, the coupling between states localized in opposite 
edges is better described using a tight-binding Hamil- 
tonian, which takes into account the coupling between 
inequivalent Dirac points. 

Electronic conductance. Armchair nanoribbons. As 
discussed above, the Dirac Hamiltonian describes ap- 
propriately the low energy band structure of armchair 
nanoribbons. Therefore, we calculate the conductance of 
the system by matching the eigenf unctions of the Dirac- 
like Hamiltonian. Given an incident electron coming 
from the left monolayer ribbon and with energy E, we 
compute the transmission coefficient to the right mono- 
layer lead. The boundary conditions on the wavefunc- 
tions determine the value of the transmission. In the 
bottom-bottom configuration (1— >1) the wavefunctions 
of the bottom layer <f>\ and <p B should be continuous at 
the beginning (x — 0) and at the end (x — L) of the 
bilayer region. For the top layer the wavefunction should 
vanish in one sublattice at x — and on the other sub- 
lattice at x = L, e.g., 4> 2 A {x = 0) = (j) B {x = L) = 0. In 
the bottom-top configuration the bottom wavefunctions 
and the top wavefunctions </>^ should be continuous 
at x — and x — L respectively. In addition, the hard- 
wall condition should be satisfied, (j> B (x = 0) = 4>\{x = 
L) = 0. Zigzag nanoribbons. To describe adequately 
the low energy properties of zigzag nanoribbons in the 
full Brillouin zone it is necessary to use a tight-binding 
Hamiltonian. We use a Green's function approach to ob- 
tain the transport properties [H [151 [05] . In this method 
the system is divided in three parts, namely, a finite-size 
bilayer section connected to the right and left monolayer 
semi-infinite leads. The Green's function of the central 
region is 

G c (e) = (e-h c -j: l -j: r )- 1 , (2) 

where He is the bilayer Hamiltonian and Ex and E# are 
the selfenergies at the ends of the bilayer region due to 
the presence of the leads. The selfenergies contain the 
information on the type of connection, i.e., 1— >1 or 1— s-2, 
of the system. In the linear regime, the conductance is 
given by 

G = 2jT(E) = 2jTr[r L g c T R g+] , (3) 

where T(E) is the transmission at the Fermi energy E 
and Tl and T B are the couplings between the bilayer 
and the left and right monolayer leads respectively. 

Results. Armchair nanoribbons. In Fig. [4] we plot the 
transmission as a function of the incident energy and the 
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FIG. 4. (Color online) Transmission as a function of energy 
and applied voltage for armchair nanoribbons in the 1— »1 and 
1— >2 configuration. The value of the interlayer hopping pa- 
rameter is 71 = O.lt. The length of the bilayer region is 
L — 35\/3et- The values of the gap edges ±Ai/2, ±A2/2 and 
± V/2, are plotted with white lines, see Fig|5] . 



applied gate voltage for an armchair nanoribbon system 
in the 1— >1 and 1— >2 configurations. The length of the 
bilayer flake is L — 35y/3a and the value of the interlayer 
hopping is 71 = t/10. The results are independent of the 
width of the ribbon, provided that the monolayer ribbons 
are metallic and that the energy of the incident electron 
is lower that the energy of the second subband. The 
transmission is obtained in the continuum approxima- 
tion, but we have checked that the results coincide with 
those obtained within the tight-binding approach. Due 
to the symmetry of the contacts, the 1— >2 configuration 
shows electron- hole and V — > — V symmetry. This is not 
the case for the 1—^1 configuration, for which the loca- 
tion of the contacts precludes those symmetries. In both 
cases the conductance is suppressed for energies in the 
gap E < |Ai|/2, for which there are no available states 
for the conductance in the bilayer region. We first discuss 
the results for the bottom-bottom configuration. In the 
energy window Ai/2 < \E\ < \V\/2 there are two prop- 
agating states in the bilayer part and the conductance 
is finite. In the range of energies 11^1/2 < l^l < A2/2 
there is only one propagating mode in the central region; 
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but in this configuration, when E and V have the same 
sign, this mode is mostly located in the opposite (top) 
layer to the leads (bottom), as can be seen in Fig. [lfb), 
and the conductance is near zero. Thus, by applying a 
gate voltage between the two layers we can tune the elec- 
tronic density in the bilayer. Changing the distribution 
of carriers from one layer to the other allows to control 
the conductance of the system by means of an external 
parameter. For energies \E\ > A 2 /2 the transmission is 
finite with antiresonances associated with interferences 
in the bilayer region due to the existence of two propa- 
gating channels. [13] These interferences are weaker for 
voltages \V\ > 71, with and overall nonzero conductance, 
because in this case the incident current from the left 
electrode is transmitted efficiently to the upper branch 
of the bilayer dispersion relation, with bottom charac- 
ter, and from there to the right (bottom) lead, with an 
almost perfect wavevector matching 20J. The weak in- 
terferences are due to the bilayer-confined states arising 
from the coupling to the top layer flake. Note the linear 
dependence of the position of the antiresonances on the 
applied voltage: the energy of the confined states in the 
top layer are displaced by the applied bias —V/2, thus 
changing the occurence of the antiresonances correspond- 
ingly. 

In the bottom-top configuration the conductance is not 
suppressed for \V\/2 < \E\ < A 2 /2 because in this case 
the incoming and outgoing electrons belong to different 
layers: the propagating mode in the bilayer has a pre- 
dominant top character (see Fig. [I]), being easily trans- 
mitted to the right electrode. For this configuration, the 
transmission at energies \E\ > A 2 /2 is generally sup- 
pressed, even though there are two propagating modes 
in the bilayer. This can be understood by noticing the 
wavevector mismatch [20] between left and right elec- 
trodes produced by the applied bias, as depicted in Fig. 
[T] Away from the gap the transmissions in the 1— >1 
and 1— >2 configurations are rather complementary; the 
antiresonances that occur in the 1—^1 configuration be- 
come resonances in the 1—^2 case. This complementarity 
of the conductance can be understood by resorting to a 
simple non-chiral model. Consider an incident carrier, 
with energy larger than the gap, coming from the left 
and therefore in the bottom sheet. When arriving at 
the bilayer central region, the incident wavefunction de- 
composes into a combination of the two eigenstates of 
the biased bilayer. The conductance through the bilayer 
region is proportional to the probability of finding an 
electron at the top (bottom) end of the central region 
for the bottom-top (bottom-bottom) configuration. As 
the total probability of finding the electron at the end of 
the bilayer region is unity, the bottom-bottom and the 
bottom-top transmissions should be the opposite. 
Zigzag nanoribbons. In Figj5]we plot the transmission as 
function of the Fermi energy and the applied gate volt- 
age for zigzag nanoribbons in the 1— >1 and 1— >2 con- 
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FIG. 5. (Color online) Transmission as a function of energy 
and applied voltage for zigzag GNRs in the 1— >1 and 1— >2 con- 
figurations. The interlayer hopping parameter is 71 = O.lt. 
The number of atoms across the nanoribbon is 16, correspond- 
ing to a width W = lla/v3- The bilayer region length is 
L — 60a. The values of the gap edges ±Ai/2 are plotted in 
white lines, see Fig. [2] 



figurations. These results have been obtained using a 
tight binding Hamiltonian and a recursive Green's func- 
tion technique [3] [T5J [IS]. The conductance of zigzag 
grapheme nanoribbons depends on the ribbon width. The 
results presented in Fig|5]correspond to a narrow ribbon, 
W = 11a/ V3, for which there is only one channel com- 
ing from the left contact for all the plotted energies. As 
for the armchair-based systems, there is a strong com- 
plementarity between the 1— >1 and the 1— >2 configura- 
tions, yielding a very different conductance as a function 
of the energy and bias for the two configurations. Other 
features, as antiresonances (resonances) in the bottom- 
bottom (bottom-top) configuration are similar to those 
occurring in armchair nanoribbons and have the same 
origin. In the upper panel of Fig. [5j corresponding to 
the 1— >1 configuration, the previous gapped region in the 
armchair case between V/2 and A 2 /2 has shrunk to a line 
of slope V/2. This is easy to understand by observing 
the zigzag bandstructure of Fig. [3j Another remarkable 
feature in Fig. [5] is the existence of a transport gap A 
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smaller than the bulk gap Ai of the gated bilaycr. As 
mentioned above, the gap A appears because of the cou- 
pling between states with the same valley polarization 
localized in different edges and moving in opposite direc- 
tions. The penetration length t of these surface states is 
rather large; for nanoribbons narrower than t this pro- 
duces a noticeable transport gap smaller than the bulk 
gap. 

In summary, we have calculated the conductance of 
bilayer graphene flakes with monolayer nanoribbon con- 
tacts with a bias voltage between layers. Depending on 
the position of the electrodes and on the applied bias, 
there is a strong variation of the conductance. Besides 
the energy gap opened by the bias, the conductance can 
be tuned by changing the spatial distribution of the car- 
riers in the bilayer region, thus allowing for the external 
control of the transport through graphene bilayer flakes. 
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